Generalized scaling theory for critical phenomena 
including essential singularity and infinite dimensionality 
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We propose a generic scaling theory for critical phenomena that includes power-law and essen- 
tial singularities in finite and infinite dimensional systems. In addition, we clarify its validity by 
analyzing the Potts model in a simple hierarchical network, where a saddle-node bifurcation of the 
renormalization-group fixed point governs the essential singularity. 



The scaling theory for the power-law singularity (PLS) 
of second-order transitions is among the most significant 
achievements in theoretical physics [H, and provides a 
comprehensive understanding of critical phenomena with 
minimal assumptions. The key concept of this theory is 
an invariance of the singular part of free energy, ga , which 
is scaled by a factor b as 



9d{t,h 



b- d g d (tb Vt ,hb v \L- l b). 
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Here, d is the dimension of space, L is the linear dimen- 
sion of the system, and t and h are deviations of model 
parameters from their critical values (such as reduced 
temperature t = T — T c and magnetic field h = H). For 
the characterization of the criticality, y t and yh are the 
most fundamental quantities, which lead to critical ex- 
ponents for the PLS of physical quantities e.g., v = 1/yt 
and j3 = (d — yh)/yt, corresponding to the correlation 
length as £ cx t~ v and the order parameter as to oc t" . 
The scaling theory is based on the self-similarity owing 
to the divergence of the correlation length. 

There is a kind of continuous transition on which stan- 
dard scaling theory does not work. The most famous ex- 
ample is the Berezinskii-Kosterhtz-ThoulessfBKT) tran- 
sition in the two-imensional XY model [2H4|, which ex- 
hibits an essential singularity (ES) around the transi- 
tion point in the disordered phase, where the correla- 
tion length diverges as £ tx exp[(t/to) -1 / 2 ] instead of the 
power-law. Below the transition temperature, £ remains 
divergent and there appears quasi-long-range order with 
zero magnetization as if the system stays in criticality. In 
this sense we hereafter call this phase the critical phase. 

A similar singularity has been recently found in infi- 
nite dimensional systems. The ES for the infinite dimen- 
sional system appears in the ordered phase; the order 
parameter behaves as to oc exp[— (i/io) -1 / 2 ] toward the 
transition point, above which the critical phase emerges. 
This is in contrast to the ordinary BKT transition, and 



is called the inverted BKT singularity Here we use 
the word "infinite dimensional" for the property that the 
typical path length L increases with system volume N 
as L oc \ogN instead of N x / d . This property is realized 
in trees, random graphs Q, hierarchical lattices hy- 
perbolic lattices pj , and small- world networks @ among 
others. Such systems have been extensively studied in the 
context of complex networks; the heterogeneous and hi- 
erarchical structure has unexpectedly revealed important 
concepts in physics Q other than the mean-field behav- 
iors for simpler infinite dimensional systems such as the 
complete graph. In particular, the ES is found rather 
often in various models on various infinite dimensional 
graphs SEMI. Therefore, the ES is considered to be 
a basic concepts of the dynamics in infinite dimensional 
systems rather than an exotic topic as the BKT transition 
in finite dimensions. Since inverted BKT singularity is 
observed in percolation U 3] and Ising models [1, H3] , 
we do not expect its mechanism to be common to that 
of the ordinary BKT transition, such as vortex-pair con- 
densation 0, y] . Universal understanding of the ES in 
infinite dimensional systems is still missing. 

In this Letter, we propose a generic scaling theory 
for critical phenomena with ES in infinite dimensional 
graphs. We also perform renormalization group (RG) 
analysis of the Potts model in a simple hierarchical net- 
work, which confirms the validity of the scaling law. 

Now, we consider a scaling formula for infinite dimen- 
sional graphs as 



(2) 



Here, 6 is a scaling factor for the total volume (mass) N 
and the correlation volume £ of the system. The length- 
based expression is obtained by replacing N — > L d , 
b — > b d , and — > £,(t) d . The volume-based expression 
can be used for the system where distance is not well 
defined, which is usual in infinite dimensional systems. 
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Equation ([2]) can be applied to various singularities by 
choosing a proper function form of A similar idea 
is seen in the finite size scaling method proposed by Kim 
fl8L fljj [it uses correlation length | (t) directly observed 
in advance]. If we assume £ is a power function as t~ x / yt , 
Eq. ((2} leads to 



ff (t,/j, A^ 1 ) = bg(tb Vt ,hb Vh ,N- 1 b), 



(3) 



which is essentially the same as Eq. ([1}. On the other 
hand, by assuming an exponential function as £(£) = 



exp[(i/i c 



we obtain 



g^h,^ 1 ) = fe- 1 ff(e- (i/to) xt b,hb yh ,N~ x b). (4) 

Note that the present formula is applicable to normal 
phases, i.e., ordered or disordered phases, but not to crit- 
ical phases where the correlation volume diverges. 

We can calculate various physical quantities by differ- 
entiating the above free energy. Hereafter, we focus on 
the case of Eq. (j4]) for the ES in the ordered phase. The 
order parameter is given by 

mit^^N' 1 ) =g h (t,h,N~ 1 ) = dg(t, h, N^/dli 

= b- {1 - yh) 'g h (e-( t/tor:ct b,hb Vh ,N^b). (5) 

First, we consider the thermodynamic limit, N = oo. 
By setting b = e^*/* - 1 "* , we have 

m = e -Ci-v*)C*/*o)— * ^ fc (! f ke^t/to)-' j ). (6) 

Since m is independent of h for h — > and independent 
of t for t — > 0, the scaling function gu(l, x, 0) should have 
asymptotic forms as 

f const, for i < 1 
g h (l,x,0) = \ , 7 

|_ x y h 1 tor x S> 1 

to reproduce 

( e -(i-2M)(V*or xt f or fay^t/to)-** < x 
to oc < . , „ . • (8) 



for he Vh ^ to "> x ' > I 



At t = 0, it reads as ma h 1 / 8 with 1/5 = y^ 1 — 1. 

Second, we consider the case of ft = with N finite, 
where 

m = e -V-v*W*°)-' t § h (l 1 0,N- 1 e<- t M~" t ). (9) 
The scaling function should be 



&(l,0,a:) 



const, for i«l 



in a similar manner as Eq. ([7]) to reproduce 

f e -(l-l/h)(t/to) -<t « f or jy-l e (t/«o) _a,t <g l 
m ^ | i\T-(l-!/h) for Ar-i e (*/*or xt > i 



(10) 
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FIG. 1: Recursive construction of a shortcut network. The 
number of generation, n, equals 1,2,3 and 4 from left to right. 
The vertical lines and solid arcs indicate backbone edges and 
shortcut edges, respectively. A periodic boundary condition 
is imposed in the vertical direction. The graph with n = 1 is 
a single node with a self-connecting edge. 



Similarly, susceptibility behaves as 

C e (2y h -l)(t/t )-*t for N -l e (t/t )-** < l 



X °C 



for JV- 1 e(*/*o) xt > 1 



(12) 



The m-th derivative of the free energy with h, which is 
proportional to £(t)( mj/,,_1 ), diverges at t = for to > 
Vfr 1 - On the other hand, the derivative with t, such as 
specific heat, never diverges. Free energy for h = and 
TV = oo is proportional to C(t)^ 1 , and the dominant term 
of its m-th derivative, t~ m ^ Xt+1 ^(t)~ 1 , goes to zero for 
t->0. 

Finally, we consider the case of t = 0. By setting b = N 
in Eq. ((5]), we have 



m = N-^- yh ^g h (0, hN y \ 1). 



The scaling function should be 



9h (0,x, 1) = 



x 



X v h 



for i<1 
for x 1 



to reproduce 



to oc 



h N 2y h -l for ft jy» h < j 

ftffc 1 - 1 for ftiV?"' > 1 



(13) 



(14) 



(15) 



where we assumed linear susceptibility for a finite size 
system. This form is the same as that of the conventional 
PLS. 

Next, we examine the validity of the present scaling 
ansatz by considering a hierarchical small- world network 
that is constructed in a recursive manner as shown in 
Fig. [TJ The graph with n generations has N = 2 n_1 
nodes and 3 x 2 n_1 — 2 edges. The degree distribution 
function is exponentially decaying as oc 2~ k / 2 . We 
note the edges represented by the vertical lines in Fig. Q] 
as backbone edges (BBEs) and arcs as shortcut edges 
(SCEs). 
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FIG. 2: (color online) Phase diagram for q = 8. The shaded 
region indicates the critical phase. The (blue) solid line and 
the (red) dashed line denote the stable and unstable fixed 
lines, respectively. The circle symbol denotes the SNB point, 

0'sn, &sn)- 



The energy function of the g-state Potts model on the 
network under magnetic field H is 

- E/k B T = E<i,j)eBBE l K + DS^o] (W, 

+ E<i,j) GS CE J5 °i°j + Ei H8 ai0 , (16) 

where 8 a b is the Kronecker's delta and Ui is a spin vari- 
able at i-th site taking one of the values, 0, 1, • • • , q — 1. 
The first summation is over BBEs and the second one 
is over SCEs. We consider only the case D = 0, but D 
becomes finite in the real space RG performed below. In 
the following we consider the case with q > 3 [22j . 

We can calculate the partition function of this sys- 
tem in a stepwise manner; decimating the spins in the 
youngest generation, which is the inverse procedure to 
grow the graph. The partial sum of the partition function 
preserves the function form by replacing the parameters 
of the energy function as 



(K n + l +D n + l S a . Q )S„. a .+(H n+l /2)(S ai a+S a . ) 



V- e {K n +D n 5, k o+J){S^ i+ S ak ^)+H n S ak0 _ £ 17 j 



where cr, and (jj are nearest spins that are older than au 
by one generation. This gives us the following recursion 
relations: 



.9n+i : 

Kn+l z 

H n +i z 
D n +\ z 

with g = 0, K 

„^< 1) - 



g n + A^/2 n+1 



(18) 
(19) 

H n +2(A^-A^) (20) 
A^-A^-2{A^-A^) (21) 

K, Hq = H, and Do = 0. Here 

e 2(K n+ D n + J) + H n + qi (22) 



e 2(K n + J) + &Hn 
„K n +D n +H n +J 



2e K n+ J e H n 



- e 
Q3 



?2 

K n + J 



<12 



(23) 
(24) 
(25) 



and q m = q — m. The quantity g n = 2~ n In C' n is regarded 
as the free energy per spin of the system with n genera- 
tions. Note that J does not change in this procedure. 

In case of no magnetic field, H = 0, both H n and D n 
remain zero and Eq. Q19p is rewritten as 



jk„(q2jk„ +2) 
qifk 2 , + 1 



where we put k n = e Kn and j = 
(FP) is obtained from k n+ i = k n 



(26) 



The fixed point 



k* as 



2qi 




(27) 



Figure [2] shows the RG fixed point and the phase bound- 
ary in k vs j space for q = 8. This system exhibits a phase 
transition from the ferromagnetic phase corresponding to 
the FPs with K = oo to the phase corresponding to the 
FPs with finite K by increasing j at fixed k. We call 
the latter 'critical phase' in the sense that the RG flow 
goes to neither K = nor K = oo but to nontrivial fixed 
points. The FP exhibits a saddle-node bifurcation (SNB) 



located at (j sn , k s 



= ( 4« 



1 



3k. 

4<ji 



q-2 

2qi 



. For k < fc s 



the phase boundary is given by a line consisting of the 
unstable FPs, each of which leads to a PLS. On the other 
hand, a transition with ES occurs at j = j sn for k > fc sn 
and its singularity is governed by the SNB point . In the 
following, we consider the phase transition for k — k sn in 
increasing j (decreasing J). 

First, we perform linear instability analysis around the 
SNB point. Equations (fT9|) - (|2~T|) arc approximated as 



(28) 



Kn + l 






D n +i 


) = M ( D n 




H n +i 


1 \H n J 







y{P\y - 


P2) 


-P2) 


(l-y)[P!(U 


' y) ' 



,(29) 



2P 2 



2P 2 (l-j/) + l 



with P3 = Pi — P-2V — 1/2. The largest eigenvalue is A+, 
and thus H n grows as e nyh with yh = lnA+. 

For H = 0, only A = 1 is the relevant eigenvalue, which 
means that the instability is marginal as expected at the 
SNB point. Equation (j2l))) is rewritten as 



in the lowest order of s n 

t - Jsn J , A — qij i ak 'l n + 1 



A{sl + Bi) 

and t, where s Tl 
and B - 



(30) 

- k 

IV r) 1 
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Equation (|30|) can be approximated to the so-called 
Kosterlitz equation Q, ds/dn = A[s 2 + Bt], which is 



(a) 



solved with s n (t) = v Bt tan 



An-JBt 



for the initial con- 



dition sq = 0. In this expression [valid until s n becomes 
0(1)], s n divergently grows as 



Sn (t) = 1/A(n {t) - n) with n (t) = ir/2AvBt. (31) 

Now let us consider the scaling behavior of the present 
model. From the aforementioned analysis, the free energy 
behaves as 

g(s = (),<;,#, TV -1 ) = 2~ n g (n (t) - n,t,He Vhn ,N~ 1 2 n ) . 

When we assume that the second argument t does not 
yield any singularity, we obtain Eq. @ by considering 
that 



, nn , , 1 , 1 /7rhi2 

= 2 , y h = ln\ + , x t = - and t = — I 



where we use 2 n °W - " = The value x t = 1/2 

must have certain universality, because it is directly de- 
rived from the simplest nonlinear RG equation, ds/dn oc 
■s 2 + Bt. This equation is a consequence of the SNB at 
the edge of the stable fixed line. We consider that this 
structure of the RG flow is the essence of the ES. In fact, 
it has been found in some systems [!, Il5l - [l7| . 

Finally we confirm the scaling ansatz by calculating 
the order parameter, m n = (rh n — q~ l )/{l — q^ 1 ) where 
m„ = % n _! + Z n -!]/dH and Z n = e 2(K n +J+D n+ H n ) + 
qi [e 2 ^ Kn+J ^+2e Hn +q 2 ] . For this aim, we additionally cal- 
culate the derivatives of g n , K n , D n and H n with respect 
to H, whose recursion equations for these quantities are 
obtained by differentiating Eqs. (|18l) - (|2"Tj) . Similarly, we 
calculate the susceptibility by using the second deriva- 
tives. 

Figure [3ja) shows the scaling plot at TV — > oo corre- 
sponding to Eqs. and (JHJ). Although we found a little 
correction to scaling that tends to disappear for H — > 0, 
a good collapse of data is obtained without any fitting 
parameter. In Fig. [3jb), we can see excellent scaling be- 
havior for t = j — j sn = corresponding to Eqs. (fT4| and 
(fT5|) . In both cases, the asymptotic form agrees with the 
prediction. 

In conclusion, we have proposed a new scaling theory 
for the ES in infinite dimensional systems, and clarified 
its validity by analyzing a simple model. We believe that 
the present scaling formulae can be applied to various 
other models, and will clarify the existence of the uni- 
versal mechanism for them. We have already confirmed 
that the scaling law holds in the bond percolation model 
on the decorated (2,2)-flower [151 [l6| (not shown here) 
and random growing network [201 ] . The finite size scal- 
ing formula included in the present theory will be useful 
in analyzing real- world data or numerical simulations, in 
which we can treat only small generations, n ~ log N. 
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FIG. 3: (color online) Scaling plot of the order parameter for 
1/N ~ (a) and for t = (b). In the latter case, we show 
the data for n=4096 and 16384 to confirm the convergence to 
the large size limit. 
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